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Experiments indicating the excitation of whispering gallery type electromagnetic modes by a vortex 
moving in an annular Josephson junction are reported. At relativistic velocities the Josephson vortex 
interacts with the modes of the superconducting stripline resonator giving rise to novel resonances 
on the current-voltage characteristic of the junction. The experimental data are in good agreement 
with analysis and numerical calculations based on the two-dimensional sine-Gordon model. 
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Whispering gallery modes are universal linear exci- 
tations of circular and annular resonators. They have 
first been observed in form of a sound wave traveling 
along the outer wall of a walkway in the circular dome 
of St. Paul's Cathedral in London and were investigated 
by Lord Rayleigh and others 0]. In the 2 meter 
wide walkway, which forms a circular gallery of 38 me- 
ter diameter about 40 meters above the ground of the 
cathedral, the whispering of a person can be transmitted 
along the wall to another person listening to the sound 
on the other side of the dome. The investigations by 
Rayleigh led to the conclusion that the whisper of a per- 
son does excite acoustic eigenmodes of the circular dome 
which can be described using high order Bessel functions. 
This acoustic phenomenon lends its name "whispering 
gallery mode" to a number of similar, mostly electro- 
magnetic excitations in circular resonators. Whispering 
gallery modes are of strong interest in micro-resonators 
used for ultra small lasers || . Most recently, circular res- 
onators with small deformations, in which chaotic whis- 
pering gallery modes were observed, attracted a lot of 
attention Q . Here we describe the experimental observa- 
tion of electromagnetic whispering gallery modes excited 
by a vortex moving in an annular Josephson junction of 
diameter less than 100 /im. 

A long Josephson junction is an intriguing nonlinear 
wave propagation medium for the experimental study 
of the interaction between linear waves and solitons ||. 
In this letter we report the excitation of whispering 
gallery type electromagnetic modes by a topological soli- 
ton (Josephson vortex) moving at relativistic velocities 
in a wide annular Josephson junction. We make use 
of the same Josephson vortex for both exciting and de- 
tecting the whispering gallery mode. These modes are 
manifested by their resonant interaction with the mov- 
ing vortex, which results in a novel fine structure on the 
current-voltage characteristic of the junction. Our ex- 
periments are consistent with the recently published the- 
ory H based on the two-dimensional sine-Gordon model. 
We present numerical calculations based on this model, 
which show good agreement with experiments. 

Electromagnetic waves in an annular Josephson junc- 




FIG. 1. Geometry of the annular Josephson junction. The 
direction of the angular velocity Q. of the vortex ($o) under 
the action of the bias current is indicated. 

tion are described by the perturbed sine-Gordon equa- 
tion (PSGE) for the superconducting phase difference (f> 
between the top and bottom superconducting electrodes 
of the junction ||. The Josephson vortex, often also 
called fluxon, corresponds to a twist over 2?r in </>. It car- 
ries a magnetic flux equal to the magnetic flux quantum 
$o = h/2e = 2.07 10" 15 Vs. Physically, this flux is in- 
duced by a vortex of the screening current flowing across 
the junction barrier. Linear excitations in this system 
are Josephson plasma waves that account for small am- 
plitude oscillations in <f>. The maximum phase velocity of 
electromagnetic waves in such a junction is the Swihart 
velocity given by cq = Xj u> p , where Xj is the Josephson 
length and lj p the plasma frequency [0. In zero exter- 
nal magnetic field the PSGE for an annular Josephson 
junction of width w < Xj can be written as 

(1 > 

where space and time are normalized by Xj and cjp 1 , 
respectively. In Eq. (0) V 2 - d 2 /dt 2 is the D'Alembert 
wave operator, sm(f> is the nonlinear term due to the 
phase-dependent Josephson current and 7 is the nor- 
malized bias current. The damping terms ad<p/dt and 
(3 V 2 d<j)/dt are inversely proportional to the quasiparticle 
resistance across the junction barrier and to the quasipar- 
ticle impedance of the electrodes, respectively. For the 
junctions of width w < Xj considered in this paper, a 
homogeneously distributed bias current 7 as in Eq. (|l|) 
is justified. In contrast, for junctions with w > Xj the 
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TABLE I. Geometrical parameters of the annular Joseph- 
son junctions used in experiment. The outer radius of every 
junction is r e = 50 /im. w = r e — n is the junction width. 



junction 


A 


B 


C 


D 


E 


n H 


47 


45 


42 


35 


30 


w [jLtm] 


3 


5 


8 


15 


20 


3 = Ti/r e 


0.94 


0.90 


0.84 


0.70 


0.60 


£ = Co/cq 


0.70 


0.78 


0.85 


0.91 


0.94 



bias current may contribute to the boundary conditions 
ofEq. (§0. 

A vortex steadily moving at a velocity u driven by 
the Lorentz force due to the bias current 7 generates 
a voltage V oc u across the Josephson junction. This 
voltage can be monitored in experiment. The radiation 
associated with the time-dependent fields described by 
Eq. (0), i.e. the magnetic field H oc |V(/>| and the electric 
field E oc d(f>/dt, can be measured either directly (for 
certain junction geometries) or through its interaction 
with the moving vortex. 

In contrast to most of the previous experiments fo- 
cusing on quasi-one-dimensional annular Josephson junc- 
tions, we investigate comparatively wide, effectively two- 
dimensional junctions. We have fabricated a set of 5 
annular Josephson junctions (A ... E) with the ratio 
5 = Ti/r e between the inner radius and the fixed 
outer radius r e = 50 /im being varied between S = 0.94 
and S — 0.60 (see Table |). The junctions are made at 
Hypres Inc. || using Nb-Al/AlO^-Nb trilayer technology 
and employ the standard biasing geometry j^] as shown 
in Fig. |l|. Due to the fabrication technology, the junction 
area is surrounded by a small passive region about 2 /im 
wide, which is omitted from Fig. [I] for clarity. In the 
passive region the top and bottom electrodes are sepa- 
rated by a 200 nm thick SiC>2 layer, which act as a small 
stripline in parallel to, but with electrical parameters dif- 
ferent from the junction itself. 

All junctions show a homogeneous bias current dis- 
tribution, inferred from the large value of the vortex-free 
critical current at zero field, which is close to the theoreti- 
cal limit. Their critical current density is j c ~ 160 Acm -2 
and the London penetration depth is Xl ~ 90 nm at 
4.2 if . The thicknesses of the top and the bottom su- 
perconducting electrode are both well in excess of Al- 
Accordingly, the characteristic parameters are estimated 
as A j « 30 /mi and v v = lo p /2tt ss 50 GHz. All pre- 
sented measurements were done at T = 4.2 K using a 
well shielded low noise measurement setup. 

We could realize single and multiple vortex states re- 
peatedly and reproducibly in any of the junctions. Vor- 
tices were trapped by applying a small bias current dur- 
ing cooling down from the normal to the superconducting 
state. Single-vortex states are identified as the lowest 
quantized voltage step observed on the current-voltage 
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FIG. 2. Experimental normalized current-voltage charac- 
teristics of single-vortex states in junctions A to E. An en- 
largement of the high voltage region of the resonances in junc- 
tions D and E is shown in the inset. 

characteristics. Also, a characteristic change of the criti- 
cal current modulation with magnetic field accompanied 
by a suppression of the critical current by a factor of 
more than 100 at zero field, as reported earlier (l(J, was 
observed when a vortex was trapped in the junction. 

In Figure || the single- vortex characteristics of the junc- 
tions A to E are shown. The current scale is normalized 
by the flux-free zero-field critical current of each junc- 
tion. The voltage of each characteristic is multiplied by 
a factor £ = co/c (see Tab. ||) calculated using the ap- 
proach by Lee et al. pT| , in order to substract the effect 
of a small passive region, cq (co) is the wave velocity in 
the junction neglecting (including) the passive region. 

A striking novel feature noticed in Fig. || is the fine 
structure on the vortex resonance which appears with 
increasing the junction width. The fine structure is most 
clearly visible for the widest junction E (see inset of 
Fig. ||) . We argue that the observed fine structure can be 
well understood as due to the interaction of the moving 
vortex with the linear whispering gallery modes Q . The 
linear modes of the annular Josephson junction resonator 
are solutions to the inhomogeneous D'Alembert equation 
in the polar coordinates (r, 9) 

{rfr r fr + ^W~ W" 1 )^ (2) 

which is found from Eq. (Q) neglecting all perturbations 
(7, ad/dt, [3V 2 d4>/dt) and approximating the nonlin- 
earity as sine/) <j) to take into account the gap in the 
plasmon excitation spectrum. In zero external magnetic 
field the boundary conditions 

^— (r = r u r e ) = (3) 

have to be fulfilled. In terms of the electromagnetic waves 
in the junction, Eq. (H) corresponds to a total internal 
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voltage, V = * a a> p / (2 jc) [^V] 
FIG. 3. Numerically calculated current-voltage character- 
istics V(y) for junctions A to E. In the inset the character- 
istics of junctions D and E are shown on an enlarged scale. 
Arrows indicate the bias points used to obtain the phase pro- 
files shown in Fig. ^[ 

reflection condition. For large angular wave numbers k 3> 
1 and wide junctions S <C 1 a solution to (0) is given by 



(f>1 ln \r,9,t) = AJ k (Lo k r)exp(ik9)ex-p(iu) k t) 



(4) 



where A is an arbitrary amplitude factor, J k is the Bcssel 
function of the first kind and oj k is the angular frequency 
associated with the mode k satisfying the boundary con- 
dition (||) at the external radius. By calculating the de- 
pendence of u>k on fc, one obtains the dispersion relation 
of whispering gallery modes in the annular junction. 

The resonance condition between the angular fre- 
quency f2 of the vortex rotation in the ring and the char- 
acteristic frequency Wfc of a whispering gallery mode can 
be expressed as 



O = LUk/k 



(5) 



In a dispersion diagram Eq. (||) determines the crossing 
points between the straight dispersion line of the vortex 
w( v ) = and the dispersion curve tj( lln ) = oj^ of the 
linear modes. At low enough damping, a vortex moving 
at the angular frequency £1 comes into resonance with a 
whispering gallery mode of wavenumber k. If the spac- 
ing in f2 between the resonances for different k is large 
enough, this effect can be observed as fine structure reso- 
nances on the single- vortex current- voltage characteristic 
of a wide annular Josephson junction. The increase of the 
excited wavenumber k with decrease of the vortex veloc- 
ity is a characteristic feature of the interaction between 
the Josephson vortex and the whispering gallery modes 
of the junction. This property is clearly manifested in 
our numerical calculations presented below. 

To confirm the above interpretation of our experi- 
mental results, we performed direct numerical simula- 
tions of Eq. (0) in polar coordinates with the bound- 
ary conditions (||) using our junction parameters (see 





FIG. 4. Phase profiles 4>{r, 6) at bias 7 equal to (a) 0.5, (b) 
0.4, (c) 0.3 on the single vortex resonance of junction E and 
at 7 = 0.6 for junction D. Plotted are lines of constant phase, 
their high density corresponds to a large gradient of phase 
and hence a large magnetic field. The position of the vortex 
is indicated by an arrow. The whispering gallery modes with 
angular wavenumber k equal to (a) 7, (b) 8, (c) 9 and (d) 9 
are observed. 

Table §). We calculated current- voltage characteristics 
V(j) = $0^(7)^/(2^) and two-dimensional phase pro- 
files 4>{r, 6, t) for various bias points using a plasma fre- 
quency of w p /27T = 52.4 GHz determined from experi- 
mental data. The damping parameter a = 0.03 was 
chosen close to its estimated experimental value at T = 
4.2 K, (3 here was set to 0. The calculated characteris- 
tics for junctions A to E are plotted in Fig. ||. Clearly, 
the fine structure on the current-voltage characteristics 
of the wide junctions D and E is very well reproduced 
in the simulation. Figure [| shows the phase profiles at 
bias points on subsequent fine structure resonances of 
junction E. Evidently, a clear whispering gallery struc- 
ture is found here. With decreasing bias, the increase of 
the angular wave number of the mode from resonance to 
resonance is observed. 

Using the resonance condition (||) and the proportion- 
ality between the angular frequency of the vortex Q and 
the voltage V measured in experiment, the fine structure 
resonances can be fitted according to the formula 



Ztt 



U 2tt k 



(6) 



where k is integer. In the limit of S <C 1 and k 3> 1 
the linear mode spectrum can be analytically approxi- 
mated as ojk = r~ 1 (k + 0.808 k 1 / 3 ) 0]. More accurately 
we have calculated the values of LUk numerically solving 
Eq. (||) with the boundary conditions (|^). The best fit 
to the experimental data is found for the plasma fre- 
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voltage, V = <u o a to / (2 k) [hV] 
FIG. 5. (a) Upper part of experimental single-vortex cur- 
rent-voltage characteristic of junction E. (b) Simulated cur- 
rent-voltage characteristic of this junction for a — 0.03, (3 = 
(open circles) and a = 0.02, /3 = 0.0012 (open squares). The 
calculated resonance voltages are indicated by vertical dotted 
lines and marked by the corresponding wavenumber k. 

quency lu p /2t: = 52.4 GHz and the wavenumber k = 7 for 
the highest voltage resonance (see Fig. ||a). This value 
of k is exactly the one found for the highest resonance 
in numerical simulations (see Fig. |]a). The voltages of 
the resonances found in simulation (Fig. |^b) are close 
to the ones calculated from the dispersion relation (dot- 
ted lines). The current-voltage characteristic of junction 
E simulated for the damping parameters a = 0.02 and 
(3 = 0.0012 is shown in Fig. ||b (open squares). Clearly, 
taking into account the quasiparticle surface losses (oc (3) 
leads to a larger differential resistance of the resonance, 
which closely resembles the experimental data in Fig. |^a. 

Considering the resonance condition Eq. ([s]) and the 
dispersion of the linear modes u>k, it can be shown that 
the density of resonances in voltage and the wave num- 
ber of the lowest excited mode increase with decreasing 
junction width w — r e — n. This fact was verified by 
numerical calculations for junction D, where the lowest 
mode number excited on the top of the resonance was 
found to be k = 9 (see Fig. ^d). For very narrow rings 
no fine structure is observed in experiment and in simu- 
lation, due to the overlapping of neighboring resonances 
in the presence of damping. 

The origin of the observed fine structure has also been 
confirmed to be due to the interaction of the vortex 
with the whispering gallery modes of the junction by in- 
vestigating its dependence on temperature, number of 
trapped vortices, and external magnetic field. At high 
temperatures, no whispering gallery resonances are ex- 
cited due to the large intrinsic damping, i.e. no fine 
structure is observed. Decreasing the temperature be- 
low 4.2 K, fine structure is observed in all samples A to E; 
also the differential resistance of the resonances decreases 
with temperature. Moreover, we found that the voltages 
of the fine structure resonances scale with the number n 



of moving vortices (or vortex/anti- vortex pairs). There- 
fore, for n > 1 the fine structure gets clearly resolved in 
voltage and also more pronounced, because several vor- 
tices coherently pump the whispering gallery mode. No 
dependence of the fine structure step voltage positions on 
small external magnetic fields was noticed. We have also 
investigated more narrow annular junctions with a wide 
idle region both experimentally and theoretically . In 
this case, the spectrum of the whispering gallery modes 
(and, thus, of the fine structure) is strongly influenced by 
the geometry and the electrical properties of the passive 
region. The fine structure recently reported in Ref. |l3| ] 
appears to be consistent with our observations. 

In summary, we have presented experimental and nu- 
merical evidences for the excitation of whispering gallery 
modes by vortices moving in wide annular junctions. 
This novel effect has been observed at sufficiently low 
damping for annular junctions in a wide range of elec- 
trical and geometrical parameters. It is very robust with 
respect to small external perturbations such as variations 
in bias current density, boundary conditions or junction 
inhomogeneitics. The resonance frequencies have been 
calculated and quantitative agreement with experimental 
data and numerical simulations better than one percent 
has been reached. Thus, the vortices appear to whisper 
(generate radiation) at frequencies between 250 and 450 
GHz in the annular whispering gallery of 100 /im diame- 
ter. 
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